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THEORETICAL  ANALYSIS  OP  A  COAXIAL  HALL  ACCELERATOR 


by 

K.Lackner 

INSTITUTE  FOR  THEORETICAL  PHYSICS 
OF  THE  INN SPRUCE  UNIVERSITY 


1 lint  roduction 


Performance  of  a  magnetohydrodynamic  accelerator  is  essen¬ 
tially  dependent  on  the  value  of  the  magnetic  interaction  parameter 


L  ...  characteristic  length 
V  ...  characteristic  velocity 


Thus,  to  achieve  efficient  interaction  while  keeping  the  length 
of  the  device  reasonably  small,  plasma  accelerators  should  ope¬ 
rate  at  high  magnetic  field  strength  and  rather  low  densities. 

In  thi3  regime,  however,  the  Hall  effect  becomes  important,  thus 
limiting  the  performance  of  devices  based  on  the  classical  Faraday 
effect . 

To  overcome  these  difficulties  several  configurations  have 
been  suggested  (see  pg.458  in  /l /  and  /2/)  in  which  the  Hall  com¬ 
ponent  of  the  plasma  currents  is  used  for  the  desired  interaction. 
Probably  the  most  interesting  one  is  the  coaxial  Hall  accelera¬ 
tor  shown  in  fig.1,  which  uses  an  applied  magnetic  field  in  the 
radial  and  an  applied  electric  field  in  the  axial  direction.  This 


d«vice  has  been  first  suggested  by  Canii  et  al.  /j>/> 

The  applied  magnetic  field  together  with  the  axial  electric 
field  gives  rise  to  Hall  currents  in  the  azimuthal  direction 
which  interact  with  the  magnetic  field  to  yield  the  desired  acce¬ 
leration  in  the  axial  direction.  The  current  component  parallel 
to  the  elctric  field  will  tend  to  cause  rotation. 


Fig.  1 

Theoretical  performance  studies  for  this  device  have  been 
carried  out  by  other  authors  /"* /,  /4/»  and  especially  by  Cann 
himself  /5/.  The  present  analysis  differs  from  the  other  ones 
in  that  the  role  of  the  magnetic  interaction  parameter  is  clearly 
pointed  out,  and  that  a  closed  analytical  solution  for  the  varia¬ 
tion  in  the  plasma  parameters  with  axial  distance  is  given. 

2)  Geometry  and  basic  equations 

The  geometry  of  the  device  is  shown  in  fig.1.  If  the  dia¬ 
meter  of  the  annular  region  is  small  or  compared  to  the  radius 
of  both  the  outer  and  the  inner  cylinder,  the  problems  can  be 
analyzed  in  cartesian  coordinates.  From  the  rotational  symmetry 


7 

S 


of  the  problem  follows  that  =  0;  by  the  same  token  also  no 
z  -  component  of  the  electric  field  will  arise.  The  assumption 


dy 


=  0  is  justified  by  the  small  ratio  of 


Summarized  these  assumptions  read: 
v  =  (u,  0,  w) 


9  9 

3  ~  3  =  ° 

y  2 


diameter  of  the  annulus 
radius 


f  -  0,  jj 


S  =  (0,  By,  0) 

2  =  (E  ,  0,  0) 

v 

Ohm's  law  is  used  in  the  form: 

j  =  +  V  X  i)  -  °eTe^  x  |  +  "  6  f2  [(i  x  $)  x  2]  (1) 

B 

(see  pg.191  of  /l /)  where  f  =  —  1  - .  (N  =  neutrals,  I  =  ions, 

I  I  1TN 

e  *  electrons). 

In  components  equ.1  reads 


>x  - : - ; — 7  [0(Ex  -  "V +  v.*.] 

’•Wi' 

j  „  - 3 - _  LuB  .  <*  t  j  1 

1+a>X.ULT,r  L  y  6  e  XJ 


e'eW 

Current  continuity  then  demands  that 
div  j  =  0:  j  =  const.  =  jQ. 

The  equations  describing  the  gasdynamic  behaviour  of  the  plasma 
-  equation  of  mass,  momentum  in  the  x  and  y  direction  and  of 
entropy  then  read: 


<0  u  »  m  =  const 


(2) 


(3) 


^*u 


dv 

dx 


-  j 


z 


-  dw  . 

t"  ■&-  i 


xBy 


^  uT 


as  £  ,  ^aV^I 

v  »-  ~  +  ■ 


.f‘ 


dx 


¥2 

J 


(/.) 

(5a) 


For  a  derivation  of  the  right-hand  side  term  in  equ.(5a)  giving 
the  Ohmic  heating  see  pg .121  of  ref.  /l/.  By  making  use  of  the 
equation  of  state  of  an  ideal  gas  (i.s.  p  = ^  RT  and  S  *  lnp  - 
0^  ln$+  SQ»  where  R,  cy  and  c^  are  referred  to  one  kg.)  equ.(5a) 
can  be  converted  into 


u( 


J_  £  a?v  _  (1 

y-1  dx  y-1  (>  dx'  ' 


+  HWif 


*2 

J 


(5b) 


3)  Operating  regimes  (local  analysis) 

The  total  input  of  electrical  energy  into  the  plasma  per 
3 

nr  and  sec  is  given  by 


j  -v-P  WpTp 

3.-E.  -  ij-f- ♦  »By  -  — jx) 


'X  X 


p  =  (1  +  MgTgwJ.'CJ  f  ) 


(6a) 


Neglecting  for  the  moment  rotation,  and  defining  a  dimensionless 
parameter  E  as 


E  = 


^Te^: 

duB 


Ve 

and  calling  — p~  =  Wt 


(7a) 

(7b) 


we  can  write  this  relation  as 


-  5  * 


VEx 


2  2 
Ou  B 
_ JL 
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1  V 

*2  2  J 
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(6b) 


Thus  the  device  will  operate  as  a  generator  -  i.e.  deliver  energy 
to  the  outer  circuit  -  if  the  above  expression  becomes  negative. 
This  is  the  case  for 


0  <  E  < 


.  •*  *• 

T 


,  2 
1  -(•<*>  T 


(8a) 


Likewise  we  may  apeak  of  a  motor  if  energy  is  fed  into  tin,  plasma, 
i.e.  if 


2  2 

E  <  0  or  E  >  (8b) 

1  +  fOT 

Note  that  motor  is  not  necessarily  equivalent  to  accelerator, 
as  the  former  definition  includes  also  the  case  in  which  the  power 
input  is  exclusively  used  for  enthalpy  increase.  This  can  be  seer, 
in  fig. 2,  where  the  oriterium  8  is  plotted  in  the  E  -  Cot  plane. 

The  definition  of  the  dimensionless  parameter  E  in  equ. 

(7a)  is  especially  suitable,  as  it  allows  for  a  formal  analogy 
to  the  parameter  E*  used  in  the  theory  of  Faraday  accelerators, 
which  is  given  by 


E 


i. 

uBz 


(7c) 


Written  in  E,  all  equations  for  the  coaxial  Hall  accelera¬ 
tor  become  identical  to  those  for  the  ideal  (i.e.  tor  »  0)  Faraday 
accelerator,  if  cot  -*  co  and  E  is  substituted  by  E*. 

This  may  be  regarded  as  an  analogy  between  the  ideal  Hail 


■&HK  rpT* 


S9MVW 


'TP«r 


-  o  - 


accelerator  (^>t  — ►  oo  )  and  the  ideal  Faraday  accelerator 
(<<n  -*  0).  Of  course  however  ion  slip  will  put  an  upper  limit 
to  the  value  of  fat  actually  achievable  (see  definition  equ.?b). 


3.2  Acceleration  and  deceleration 

To  determine  the  sign  of  the  change  of  the  flow  parameters 
with  axial  distance  x  it  is  necessary  to  bring  equ.(3)  and  (5) 
into  canonic  form,  where  use  has  to  be  made  of  the  algebraic  re¬ 
lations  (l)  and  (2).  (Equ,(4)  can  be  separated  from  (3)  and  (5)). 
As  a  result  one  obtains 


1 

u  dx 

mg  M2  ( 

*  2  *  2 
dB  M  -1 

o 

[  (e-i ) 

-  (r-D  L< 

&  T 

(9a) 

1  dp 
p  dx 

_  i 

2  2 
ob:  m  -i 

o 

[(E-1) 

-  (r-i)M2j 

<0  T 

(9b) 

Thus  to  yield  acceleration,  the  expression  in  braces  at 

the  right  side  of  (9a)  has  to  be  larger  than  zero  in  the  supersor-1  * 

and  smaller  than  zero  in  the  subsonic  case.  Concentrating  on  the 
2 

former  (M  >>  1 )  we  obtain  that  for  acceleration 


2  2 
CO  T 

2  2  , 
to  1  +1 


(£x=i 

'2y-° 


3  4 


(2y-2)2  (y-l)u>V 


)<  s< 


C  i 

<0  T 


2  2 

U>  T  +  i 


(22=1  + 

' 2  y-2  + 


f 


( 2  y— 2 )  * 


(y-l)w2t2 


(10a) 


has  to  be  valid.  For  wr  -*  oo  we  obtain 


1  <  E*  <  -*•? 

y-1 


-  7  - 


the  corresponding  relation  for  the  Faraday  accelerator.  Further¬ 
more  we  can  see  from  (10a)  that  for  values  of 

6>t  <  4y(y-i ) '  ( i ot) 

acceleration  cannot  be  achieved  by  any  combination  of  the  other 
parameters  (in  the  supersonic  case).  Equ.9a  may  also  be  used  to 
determine  the  value  of  E  yielding  maximum  acceleration,  if  the 
other  parameters  are  kept  constant.  As  a  result  we  get 


2  2  „  , 

v  .  i£-L _ n 

^AX  "  i+t it2  ^*2' 


(10c) 


Both  the  boundaries  of  the  accelerator  regime,  as  well  as  Ej^x 
are  shown  in  fig. 2  for  y  =  1 . 67  *  Only  for  large  values  of  cot  the 
lower  boundary  of  the  accelerator  regime  coincides  with  the  upper 
one  for  generator  action.  This  also  follows  from  the  formal  ana¬ 
logy  to  the  ideal  Faraday  type  device. 


3. 3  Pressure  variation 

Equation  (9b)  may  be  used  to  del ive  relation  similar  to 
(10a  -  c)  for  the  variation  in  pressure.  For  supersonic  flow, 
pressure  will  decrease  downstream,  if 

iJxL_  -1/  1  i  ;  1  _  r~~y 

Jt2+1  2(y-l)M2  '  4(H)V  A>S2  (y-l)M2w2t2 

<££—  (2(,Y-..i.).M2_ti,  /  t  . 1  '  ~  ~  „  1„„  _  .  1  , 

Ov2+ 1  2(y-l)M2  »  4(y-l)2M4  fc2?2  (y-l)M2w2T2 


E< 
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Otherwise  the  pressure  will  increase.  Contrary  to  the  criterium 
for  velocity  increase,  the  bounds  for  E  given  by  relation  (11a) 
depend  explicitely  on  the  Mach  number.  Maximum  pressure-  decr;:s>. 
is  achieved,  when  E  is  chosen  to  satisfy 


2  2 

E  -  fV?  + 

1  +  Cd  t 


2(y-1 )M 


2  1 


3.4  Variation  of  Mach  number 
Making  use  of 

M2  _  m 

-  a2  "  yP  ~  yp 
which  differentiated  yields 


dM  _  1  /du  dp  v 
M  =  2  ^  u  ‘  p  ' 


(11'' ' 


equations  (9a)  and  (9b)  can  be  combined  to  give 

.2 


1  dM  _  mg  _M 
M  dx  2 OB 

o 


2  2 

(  (y-1  )(E-1 )  -  (y-l)(yM2+l)  [  (E-l)2+  -§-*■]{(  9c) 
i  M  -1  v  L  U>  T  ) 


In  supersonic  flow,  increase  in  Mach  number  thus  occurs  only,  if 
E  satisfies 


c^t2  /2yM2(y-1  )+3v 
(it2+1  2(y-l)(yM2+1 


-  -V — 

)  ’  4(r-' 


e<  -£L-  f.Ay-'hlr- 1  £ 

wS+ 1  2(y-1 )(yK  +1 )  ’  , 


iltl.ll _  .  2¥li,zzllt2r _ )  <r 

)2(yM2+l)2  w2t2(y-1 )(yM2+1 ) 

(12a) 


(ill  )2  _  yM  ( y-1  )+2y  \ 

toT^+l  2(  y-1 ) ( yK^+1 )  1  4( y-1 )2 (yM2+1 )“  W t2 ( y-1 ) ( yM2+1 ) 


i 

i 


1  ?***>., 


-  Q  ~ 


In  fig. 3  the  regimes  in  which  increase  m  Mach  number  occurs  are 
shown  for  five  values  of  M  (full  lines).  The  smalles*  v^iue  used 
is  slightly  above  1 ,  as  M  -  1  for  itself  is  a  singular  case  (see 
equ.9c) . 

From  (12a)  then  follows  that  for  each  Mach  number  a  mini¬ 
mum  value  of  the  Hal]  parameter  o t  exists,  below  which  increase 
in  Mach  number  can  never  be  achieved  (see  fig./'.  This  value  of 
is  given  by 


h)T  = 


r  +  1 


((y-D(yM  +1)  .  1  \ 

v  y  +  i  ' 


(12b) 


Finally  one  can  find  that  for  eacu  value  of  the  Mach  number  M, 
maximum  increase  of  M  will  be  obtained  for 


2  2 
1+uit 


2(y-l)(yM  +1) 


(12c) 


For  M  slightly  above  1  and  for  M  =  2  the  corresponding  curves  are 
also  given  in  fig. 3  (dashed  line). 


3.5  Optimum  power  efficiency 

The  usual  optimization  criterium  for  propulsion  applications 


is  maximization  of  the  ratio 


sower  input  via  Lorentz  force 
total  power  input 
E  -  1 


.  In  our 


parameter  this  ratio  fs  becomes  * - - -  . 

E(E(u4-y)-1) 


It  thus  has  a  maximum  for 


E  =  E 


'^17 


Also  this  relation  is  represented  graphically  in  fig. 2, 
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4)  Total  performance 

To  study  the  overall  behaviour  of  a  coaxial  Hall  accelera¬ 
tor,  the  governing  equations  (2)  -  (5)  have  to  be  integrated  with 
respec:  to  x  over  the  whole  accelerator  length.  For  this  purpose 

we  introduce  the  following  dimensionless  variables: 

2 

2  U  mUo 

G  =  u/u  ,  w  =  w/wq,  p  =  p/pQ,  m  =  e0«0»  M0  “  *T~  =  -ZT  (I4a-e) 

a  '  *  o 

2  0 
3 

as  =  dx  (i4f) 


The-  introduction  of  the  interaction  parameter  S  as  independent 
variable  makes  the  following  analysis  applicable  also  to  the  case 
of  variable  electrical  and  cf  varying  magnetic  field  strength 
without  further  modification  (of  course  always  under  the  condition 
that  the  assumptions  concerning  the  geometry  remain  valid).  On 
the  basis  of  a  kinetic  argument  given  by  Cann  /5 /  (his  quantity 
is  the  reciprocal  of  the  parameter  E  used  in  our  analysis)  ,  it 
can  be  shown  that  E  will  remain  constant  over  the  accelerator 
length;  provided  the  degree  of  ionization  is  constant  too.  To 


see  this  we  insert  the  relations: 

2 


“ —  Tg  (see  /l /  pg  .1  90) 
6 


(15a) 


j  =  const  (since  div.j  =  0) 


(15b) 


_ 


(15c) 


n_  n 

a  =  —  =  —  (assuming  quasineutrality)  (I5d) 
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nm  u  e  m  (m  being  the  atomic  mass) 

B  & 


(I5e) 


into  equation  (7a)  »  obtaining: 


B  = 


<JuB 


m  w  t  nm  j  j  m 
e  e  e  au3t  °x  a 


n  e  t  B  m 
e  e  y 


mae 


const 


The  governing  equations  then  read 


*  0 


3  at  +  rp  at  -  r(r-i)M2  [(gj)2  +  (e-i)2J  3 


(l6a) 

(16b) 


Jv  * 

dS  Wx 


u 


(16c) 


After  some  manipulations  equation  (l6a)  and  (l6b)  can  be  combined 
into  the  second  order  equation 

(17) 


52{J" 


au  1  ^  +  buu 1  2  =  0 


where 


a  = 


v  +  1 


(r-')[(^)2t(E-i)2]  -  ( e— i ) 


y(E  -  1) _ 

(y-i)[(~)2  +  (e-i)2J  -  (E-i) 


By  the  transformation 


u  = 

du 

dS 


X 

(18a) 

=  y 

(18b) 

equation  (17)  can  be  transformed  into  the  first-order  equation 


dy  a  2  b 

d?  =  7y  - 


(19) 
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which  is  of  the  Bernoulli  type  and  can  be  tranformed  into 
linear  one  by  the  substitution 


y 


(20a) 


yielding 


with  the  solution 


z>  -  -  z  =  - 

x  2 

x 


z  =  Cxb  +  ~~  1 
b+1  x 


After  re-transformation  and  after  carrying  out  the  integrate 
of  (l8b)  we  obtain  the  implicite  solution 


*  i„5  .  ^l(s-s0) 


The  constants  c'  and  SQ  have  to  be  determined  by  the  initial 


conditions 


S  =  0:  u  =  1 


p  =  1 


One  then  obtains 


:'(ub+1  -  1)  + 


ln„  .  ktl  s 

a 


(24a) 


P  -  M=[(r  ♦  +^U)  +  U1 


where 


(24b) 


M^(y+1)  b+1 

Equation  (l6c)  can  now  be  integrated  too,  to 


give 
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’  'JrtbTa  °'<ub+2  -  ’> +  bfr  <u  - ')]  <24°> 

where  use  has  been  made  of  £f(u)dx  =  ^ f(u)~-,  and  the  relation 

(18b)  was  substituted  for  u'.  Similarly  we  obtain  the  total  power 
input  via  Lorentz  force  as 

h  ■  ?ouo  [(E-')u2CS  -  f0u l  (E-1)^  ub+2  +  »2(  (244) 

These  results  are  illustrated  on  fig. 5-8.  Pig. 5  is  a  plot  of 
velocity  versus  interaction  parameter  for  =  6.0,  y  =  1,67 
and  for  E  =  2.0  (case  _)  ,  1  .69  (case  II),  1 . 1 6  (case  III).  The 
last  two  values  correspond  to  E^AX  and  E^pp  (equ.lOc  and  13). 

As  pointed  out  already  by  Cann  / 5/  velocity  is  increasing  without 
limiting  value.  The  dashed  line  gives  the  analogous  curve  for 
E  =  2.0  and  cjx  =  3-0  (case  IV). 

Fig. 6  shows  the  corresponding  results  for  the  variation 
in  pressure,  again  for  the  above  four  cases. 

Pig. 7  illustrates  the  influence  of  the  entrance  Mach  num¬ 
ber  on  velocity  (full  line)  and  pressure  variation  (dashed  line). 
Shown  are  the  cases  M  =1.2  (curve  I)  and  M  =3.0  (curve  II ); 

E  =  2.0  and  =  6.0  in  both  cases. 

Fig. 8  is  a  plot  of  the  interaction  parameter  S  necessary 
for  a  tenfold  increase  in  velocity  versus  Hall  parameter.  Curve  I 
corresponds  to  E  =  2.0,  curve  II  to  the  case  EyAX  (i.e.  for  each 
40*  the  value  E  satisfying  relation  (l  0<^  is  taken),  curve  III  to 
the  case  E0pr. 
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